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Abstract 

We study invariant manifolds of measure-valued solutions of the 
partial differential equation for geodesic flow of a pressureless fluid. 
These solutions describe interaction dynamics on lower-dimensional 
support sets; for example, curves, or filaments, of momentum in the 
plane. The 2+1 solutions we study are planar generalizations of the 
1+1 peakon solutions of Camassa & Holm [1993] for shallow water 
solitons. As an example, we study the canonical Hamiltonian inter- 
action dynamics of N rotating concentric circles of peakons, whose 
solution manifold is 2N— dimensional. Thus, the problem is reduced 
from infinite dimensions to a finite-dimensional, canonical, invariant 
manifold. The existence of this reduced solution manifold and many 
of its properties may be understood, by noticing that it is also the mo- 
mentum map for the action of diffeomorphisms on the space of curves 
in the plane. We show both analytical and numerical results. 

PACS numbers: 47.35, ll.10.Ef, 45. 20. J] 

Keywords: Shallow water flows, Hamiltonian systems, Momentum filaments. 



Contents 

ll Introduction and 



overview 



2 Measure-valued momentum maps and solutions of geodesic flow in n dimensions! 



3 Lie-Poisson bracket for rotating concentric circles of peakons! 11 



Azimuthallv symmetric peakons! 13 

4.1 Derivation of equations! 13 



4.2 Solution propertied 16 



Numerical Results for Radial Peakons! 17 

5.1 Radial peakon collisions! 17 

5.2 Bouncing off the center I 20 



6 Numerical Results for Rotating Peakon Circles! 21 

7 Conclusions! 26 



D. D. Holm, V. Putkaradze & S. Stechmann Rotating Circular Peakons 3 



8 Appendix: Extension of one-dimensional peakonsl 28 

1 Introduction and overview 
Geodesic flow in n dimensions 

As first shown in Arnold [1966] pQ, Euler's equations of ideal fluid dynamics 
represent geodesic motion on the volume-preserving diffeomorphisms with 
respect to the the L? norm of the velocity. More generally, a time-dependent 
smooth map g(t) is a geodesic on the diffeomorphisms with respect to a 
kinetic energy norm KE = |||u|| 2 , provided its velocity, the right-invariant 
tangent vector u = gg~ l (t), satisfies the vector Euler-Poincare equation, 

d t m = — u ■ Vm — Vu T • m — mdivu = — ad *m . (1) 

Here ad* is the adjoint with respect to L 2 pairing ( • , • ) : q* x q — > R of the 
ad- action (commutator) of vector fields u, w G g. That is, 

( ad * m , w ) = — ( m , ad u w ) = — ( m , [u, w] ) . (2) 

The momentum vector m G Q* is defined as the variational derivative of 
kinetic energy with respect to velocity, 

5{KE) = (m,5u) & m = . (3) 

This defining relation for momentum closes the Euler-Poincare equation (JTJ) 
for geodesic motion with respect to the kinetic energy metric KE = i||u|| 2 . 
For more details, extensions and applications of the Euler-Poincare equation 
to both compressible and incompressible fluid and plasma dynamics, see 
Holm, Marsden and Ratiu [1998] [12]. 

Geodesic flow with H 1 velocities in two dimensions 

In this paper, we consider the solution behavior of the Euler-Poincare equa- 
tion (£[]) when the momentum vector is related to the velocity by the two- 
dimensional Helmholtz operation, 



m = u — Au , 



(4) 
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where A denotes the Laplacian operator. This Helmholtz relation arises 
when the kinetic energy is given by the H 1 norm of the velocity, 



KE= -Hull^! = iy |u| 2 + \Wu\ 2 dxdy. (5) 

The H l kinetic energy norm (jSJ) is an approximation of the Lagrangian in 
Hamilton's principle for columnar motion of shallow water over a flat bottom, 
when potential energy is negligible (the zero linear dispersion limit) and the 
kinetic energy of vertical motion is approximated by the second term in 
the integral, [3]. In this approximation, the physical meaning of the quantity 
m = 5(KE) / 5u in the Helmholtz relation (j3J) is the momentum of the shallow 
water flow, while u is its velocity in two dimensions. See Kruse et al. [2001] 
j3] for details of the derivation of the geodesic equation for approximating 
2D shallow water dynamics in this limit. 

Problem statement: Geodesic flow on H l in cylindrical symmetry 

The present work studies azimuthally symmetric solutions of the Euler- 
Poincare equation (£[]) in polar coordinates (r, <fi) , 

m = m r (r, t) f + m^(r, t) cf) and u = u r (r, t) r + u<j>(r, t) <j> . (6) 

In the standard basis for cylindrical coordinates, momentum one-forms and 
velocity vector fields are expressed as, 

m ■ dx = m r dr + rm^ d(f) and u ■ V = u r d r H — -8$ . (7) 

Solutions © satisfy coupled partial differential equations whose radial and 
azimuthal components are, respectively, 

8m / Vb,k \ 

= --8 r (rm r u r ) -m r 8 r u r - (rm$)8 r [^j , (8) 

^ = -\d r {r*m,u r ) . (9) 

In these coupled equations, nonzero rotation generates radial velocity u r , 
which influences the azimuthal motion. Without rotation, = 0, and the 
solution becomes purely radial. The system of equations (jSJ) and (jHJ) for 
geodesic motion conserves the H 1 kinetic energy norm, 

#E(K],M) = ^ J [ul + ul r +^ + u 2 r + ul r + ^Ydr = ^\\u\\ 2 H1 . (10) 
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The corresponding momenta satisfy the dual relations, 

8KE . 8KE 



and m r 



The Helmholtz relation (J3J), between m$ and u^, for example, becomes 



(11) 



II,: 



d 2 u<f ) 1 du^ 1 



<9r 5 



r <9r 



-d r rd r 



u„ 



'121 



This relation between velocity and momentum defines the symmetric invert- 
ible Helmholtz operator in cylindrical geometry with finite boundary condi- 
tions at r = and — ■> oo. The velocity w(r, £) is obtained from the momentum 
m(r, t) by the convolution u(r) = G *m = J°° G(r, £)m(£)£d£ (an extra fac- 
tor £ arises in cylindrical geometry) with the Green's function G(r, £). The 
Green's function for the radial Helmholtz operator in (fT2*|) is given by, 



G(r,0 



7i(0^i(r) for £ < r 
Ji(r)iTi(0 for r<£ 



(13) 



where Ji and iTj are modified Bessel's functions. This Green's function will 
play a significant role in what follows. 

Equations (JHJ) and Q may be written equivalently as a Hamiltonian 
system, by Legendre transforming the kinetic energy KE([u r ], [u^]) to the 
Hamiltonian h([m r ], [rm^]), after which the equations take the form, 



0_ 

Of 



where Sh/Sm r 



{m r , h} 
{rm^, h} 



V 



5h/5m r 
8hj8{rm^) 



-V 



u r 



(14) 



u r , 8hj8{rm^ 
V 



u^/r and the Hamiltonian operator D, 



r l d r rm T + m r d r rm^d r 
r~ l d T r 2 m t j ) 



(15) 



is the matrix T> = {(m r ,rm^) , (m r ,rm0)}, which defines the Lie-Poisson 
bracket for geodesic motion in cylindrical geometry. Note that D is skew- 
symmetric with respect to the L 2 pairing with cylindrical radial measure. 
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2 Measure-valued momentum maps and so- 
lutions of geodesic flow in n dimensions 

Measure-valued solution ansatz 

Based on the peakon solutions for the Camassa-Holm equation [2] and its 
generalizations to include the other traveling- wave pulson shapes jH], Holm 
& Staley jH] introduced the following measure- valued ansatz for the solutions 
of the vector EP equation (0) , 

N 

m(x, t) = V / P a (s, t) 6( x - Q a (s, t))ds, m G M n , s G R k , (16) 



.*) = £ [P a (s,t)5{x-Q a (s,t))ds, 

a=l ^ s 



where the dimensions satisfy k < n. The fluid velocity corresponding to the 
momentum solution ansatz ()16j) is given by 

- r 

u(x,t) = G*m = / P b (s',t)G(x,Q\s',t))ds' , ueR", (17) 
6=1 -^ s ' 

where G(x, y) is the Green's function for the Helmholtz operator in n di- 
mensions. These solutions are vector-valued functions whose momenta are 
supported in M n on a set of N surfaces (or curves) of codimension {n — k) 
for s G M fc with k < n. In three dimensions, for example, they may be 
supported on sets of points (vector peakons, k — 0), quasi one-dimensional 
filaments (strings, k = 1), or quasi two-dimensional surfaces (sheets, k = 2). 
Substitution of the solution ansatz (j!6)l into the EP equation (0) implies the 
following integro-partial-differential equations (IPDEs) for the evolution of 
such strings, or sheets, 

B N r 

-Q a (M) =J2 / P\s',t)G(Q a (s,t),Q b (s',t))ds', (18) 

6=1 

a N r a 

-P°(M) = |,(HM)-PV,i))^G(Q fl (M),QV,t))^. 



Importantly for the interpretation of these solutions given later in Holm 
and Marsden [2003] [TT], the independent variables s G M. k turn out to be 
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Lagrangian coordinates. When evaluated along the curve x = Q a (s,t), the 
fluid velocity (fTTjl satisfies, 

ucm) =y; / p v, t) g{ q a (s, t ), q v, t) )d S > = dQa ^ t] . 

(19) 

Consequently, the lower- dimensional support sets (defined on x = Q a (s,t) 
and parameterized by coordinates s G M fc ) move with the fluid velocity. 
Moreover, equations (fTSjl for the evolution of these support sets are canonical 
Hamiltonian equations, 

1q"( s t)- 6 -^- -P a (s t)--— (20) 
The corresponding Hamiltonian function Hn : (M n x IR' 1 )®^ — > R is, 




This is the Hamiltonian for canonical geodesic motion on the cotangent bun- 
dle of a set of curves Q a (s,t), a = 1, . . ., N, with respect to the metric 
given by G. 

Hamiltonian geodesic dynamics under the measure- valued vector solution 
ansatz for the EP equation JI} was investigated numerically in Holm & 
Staley [2003b] ^U]- We refer to that paper for more details of the solution 
dynamics. 

The non-locality of the dynamics of the measured-valued solutions, de- 
scribed in the literature, makes analytical progress difficult. In this paper, we 
derive measured-valued solutions that are rotationally symmetric. For this 
circular symmetry, the non-locality integrates out and the motion reduces 
to a set of ordinary differential equations. Most of the paper is devoted to 
these circular solutions, which we call rotating peakons. The set of solu- 
tions obeying translational symmetry, but having two velocity components, 
is studied in the appendix. For a solution of N planar peakons, there are 
2N degrees of freedom, with 2N positions and 2N canonically conjugate mo- 
menta. Hence, evolution of the N planar peakon solution is governed by a 
set of 4iV nonlocal partial differential equations. These reduce to ordinary 
differential equations in the presence of either rotational, or translational 
symmetry. We also demonstrate that these solutions emerge from any initial 
condition with this symmetry in the plane. 
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Potential applications of measure-valued solutions of geodesic flow 

One of the potential applications of the two-dimensional version of this prob- 
lem involves the internal waves on the interface between two layers of different 
density in the ocean. Fig. ^ shows a striking agreement between two inter- 
nal wave trains propagating at the interface of different density levels in the 
South China Sea, and the solution appearing in the simulations of the EP 
equation in two dimensions. Inspired by this figure, we shall construct a 
theory of propagating one-dimensional momentum filaments in two dimen- 
sions. For other work on the 2D CH equation in the context of shallow water 
waves, see Kruse, et al. (2001) |3|. 

Another potential application of the two-dimensional version of this prob- 
lem occurs in image processing for computational anatomy, e.g., brain map- 
ping from PET scans. For this application, one envisions the geodesic mo- 
tion as an optimization problem whose solution maps one measured two- 
dimensional PET scan to another, by interpolation in three dimensions along 
a geodesic path between them in the space of diffeomorphisms. In this situ- 
ation, the measure- valued solutions of geodesic flow studied here correspond 
to "cartoon" outlines of PET scan images. The geodesic "evolution" in the 
space between them provides a three dimensional image that is optimal for 
the chosen norm. For a review of this imaging approach, which is called "tem- 
plate matching" in computational anatomy, see Miller and Younes [2002] ^3] • 

Peakon Momentum Map J : T*S — > Q* in n dimensions 

Holm and Marsden [2003] explained an important component of the 
general theory underlying the remarkable reduced solutions of the vector 
EP equation In particular, Holm and Marsden [2003] showed that 
the solution ansatz (fTBj) for the momentum vector in the EP equation (0) 
introduced in Holm and Staley [2003] fU] defines a momentum map for the 
action of diffeomorphisms on the support sets S of the Dirac delta functions. 
These support sets are points on the real line for the CH shallow water 
equation in one dimension. They are points, curves, or surfaces in MJ 1 for the 
vector EP equation in n— dimensions. 

Momentum map definition. Let a group G act on a manifold S, and 
lift the action of G to the cotangent bundle T*S. A momentum map J 
is a Poisson map from T*S to g*, the dual of the Lie algebra of G. (A 
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Figure 1: Simulation of the full EP equation (JTJ), courtesy of Martin Staley 
(top). Internal waves in the South China Sea (bottom). 
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map is Poisson, provided it is coadjoint equivariant. In particular, J maps 
the canonical Poisson bracket on the image space T*S into the Lie- Poisson 
bracket on the target space g*.) In symbols, this is 



J: (P,Q) ET*S — >mGg*, 

J: {/,/i} can (P,Q)^{/,/i} iP (m) = (m. 
where (•,•): 0* x g — > R . 



Sf Sh 
<5m ' <5m 



The n— dimensional peakon momentum solution ansatz J (for any Hamil- 
tonian) is given by Holm and Staley [2003] ^U] as the superposition formula 
in (US!), 



A' 



J : m(x, t) = ^2 f P a (s, t)5(x- Q a (s, t))ds, meR n ,seR k . (22) 

o=l •** 

By direct substitution using the canonical Q, P Poisson brackets, one com- 
putes the Poisson property of the map J in n Cartesian dimensions. Namely, 

{ m! (x), mj (y)} c jP,Q) = - (— m! (x)+m 3 (x)- )^(x-y), (23) 

in the sense of distributions integrated against a pair of smooth functions of x 
and y. This expression defines the Lie- Poisson bracket {• , •} i - p (m) defined on 
the dual Lie algebra g*, restricted to momentum filaments supported on the 
N curves x = Q a (s, t), where a = 1, 2, . . . , N. Its calculation demonstrates 
the following. 

Theorem 2.1 (Holm and Marsden [2003]) The momentum solution ansatz 
Mb}) , for measure-valued solutions of the vector EP equation (QJ), is a mo- 
mentum map. 



The Poisson property of the momentum map J in (J22j) is, of course, indepen- 
dent of the choice of Hamiltonian. This independence explains, for example, 
why the map extends from peakons of a particular shape, to the pulsons of 
any shape studied in Fringer and Holm [2001] [SJ. The solution ansatz (fTo^ 
now rewritten as the momentum map J in (122(1 is also a Lagrange-to-Euler 
map, because the momentum is supported on filaments that move with the 
fluid velocity. Hence, the motion governed by the vector EP equation ([T|) 
occurs by the action of the diffeomorphisms in G on the support set of the 
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fluid momentum, whose position and canonical momentum are defined on 
the cotangent bundle T*S of the space of curves S. This observation informs 
the study of geodesic motion governed by equation (JTJ. For complete details 
and definitions, see Holm and Marsden [2003] |llj . 

Peakon momentum map J : T*S — > Q* on a Riemannian manifold 

The goal of the present work is to characterize the measure- valued solutions 
of the vector EP equation (JTJ, by using the momentum map J in (J22j) when 
S is the space of concentric circles in the plane. The motion and interac- 
tions of these measure- valued solutions may be either purely radial (circles of 
peakons), or they may also have an azimuthal component (rotating circles of 
peakons). To accomplish this goal, we employ a result from Holm and Sta- 
ley [2003b] [H]] that on a Riemannian manifold M with metric determinant 
det g(x), the measure- valued momentum ansatz (|16|) becomes 

- f 6(x — O a (s t) ) 
m(x, t) = V / P a (s, t) -i JLl' ds, m E M, s E R k . (24) 

This solution ansatz is also a momentum map, as shown in Holm and Mars- 
den [2003] [TJ. On a Riemannian manifold, the corresponding Lie-Poisson 
bracket for the momentum on its support set becomes 

{m i( x) , m,(y)} = - (m,.(x) A + _L_ A_ ^m^)) ^£ .(25) 

For example, in cylindrical symmetry one has y/det g = r and the vector m 
depends only on the radial coordinate r. For solutions with these symme- 
tries, the Lagrangian label coordinate s is unnecessary, as we shall see in the 
cylindrical case, and the equations for Q a and P a will reduce to ordinary 
differential equations in time. 

3 Lie-Poisson bracket for rotating concentric 
circles of peakons 

Azimuthal relabeling symmetry for rotating circular peakons 

We shall consider the dynamics of circles of peakons, whose motion may 
have both radial and azimuthal components. These are rotating circular 
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peakons. Suppose one were to mark a Lagrangian point on the a— th circle, 
a = 1, . . . , N. Then the change in its azimuthal angle 4> a {t) could be measured 
as it moved with the azimuthal fluid velocity along the the a— th circle 
as its radius r = q a (t) evolved. Translations in the Lagrangian azimuthal 
coordinate would shift the mark, but this shift of a Lagrangian label would 
have no effect on the Eulerian velocity dynamics of the system. Such a 
Lagrangian relabeling would be a symmetry for any Hamiltonian depending 
only on Eulerian velocity. Thus, the azimuthal relabeling would result in 
the conservation of its canonically conjugate angular momentum M a , which 
generates the rotation corresponding to the relabeling symmetry of the a— th 
circle. The a— th circle would be characterized in phase space by its radius 
r = q a (t), and its canonically conjugate radial momentum, denoted as p a . 
The rotational degree of freedom of the a— th circle would be represented by 
its conserved angular momentum M a and its ignorable canonical azimuthal 
angle <p a . The only nonzero canonical Poisson brackets among these variables 
are, 

{q a ,Pb}can = 8ab and {(fi a , M b } can = 5 ab . (26) 
Momentum map for rotating circular peakons 

In terms of their AN canonical phase space variables (q a ,p a , 4> a , M a ), with a = 
1,2, . . ., N, the superposition formula (j22j) for N rotating circular peakons 
may be expressed as, 

^(.^ElrflH^*)^. (27, 

We shall first verify that this formula is a momentum map, and then in 
section 0] we shall derive it, by requiring it to be a valid solution ansatz 
for the geodesic EP equation (0) in polar coordinates. As a consequence, 
the motion governed by the system of partial differential equations (jSJ) and 
for geodesic motion in the plane with azimuthal symmetry has a finite 
dimensional invariant manifold in the 2iV— dimensional canonical phase space 
( 9a j Pa) f° r each choice of the N angular momentum values M a , with a = 
1, 2, . . . , N. Later, we shall also examine numerical studies of these solutions 
when the kinetic energy is chosen to be the H 1 norm of the azimuthally 
symmetric fluid velocity. 



D. D. Holm, V. Putkaradze & S. Stechmann Rotating Circular Peakons 13 



By direct substitution using the canonical Poisson brackets in (|2fj|) . one 
computes the Poisson property of the map J in (|27jl. Namely, 



{m r (r),m r (r')} can (p,q) = - (~m r (r) + m r (r)^j — 

{m r (r),r / m^(r / )} can (p,q) = -rm^(r)— , 

r ^ / ai < \ 1 d 2 i \ K r - r ') 

{rm^{r),m r (r )) (p,q) = - -— r m^r) 



r') 



IT C^T" f 

(rm^(r),rm (/) (r')} can (p,q) = 0. (28) 

These equalities are written in the sense of distributions integrated against a 
pair of smooth functions of r and r' . They demonstrate the Poisson property 
of the map J in (j2*7|) . which is also the solution ansatz for the rotating circular 
peakons. They also express the Lie-Poisson bracket {■ , •} iJ ,(m n rm^) for 
momentum filaments defined on the dual Lie algebra g* and restricted to the 
support set of these solutions. Hence, we have demonstrated the following: 

Proposition 3.1 The map J in p?7|) is a momentum map. 



On comparing the formulas in (j28j) with the Hamiltonian operator T> 
for the continuous solutions in (|15[). one sees that the momentum map (|27j) 
essentially restricts the Lie-Poisson bracket with Hamiltonian operator T> to 
its support set. Next, we shall derive the momentum map (J27)) by requiring 
it to be a valid solution ansatz for the geodesic EP equation (JTJ in polar 
coordinates. 



4 Azimuthally symmetric peakons 
4.1 Derivation of equations 

We seek azimuthally symmetric solutions of the geodesic EP equation (JTJ) in 
polar coordinates (?", 0), for which 

m = m r (r, t)r + m^r, t)(p = (m r (r, t),m^,(r, t)) . (29) 

We shall derive the momentum map (J27j) and the canonical Hamiltonian 
equations for its parameters (q a ,p a , M a ) by assuming solutions in the form, 

m(r,t) = jh (piW + Vi®*) *< r (30) 

8=1 
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These solutions represent concentric cylindrical momentum filaments which 
are rotating around the origin. The corresponding velocity components are 
obtained from 

(u r (r, t), U(f> (r, t)) = [ r'G(r, r') (m r (r', t), m^r', t)) dr' , (31) 



where G(r, r') = G(r', r) is the (symmetric) Green's function for the radial 
Helmholtz operator given in formula (fT3*j) . Hence, the fluid velocity corre- 
sponding to the solution ansatz ()30|) assumes the form, 



A 1 



U M) = E (pM* + <&(*)) 



(32) 



with Green's function G(r,qj(t)) as in formula (|13J). In addition, the kinetic 
energy of the system is given by 

If 1 N 

KE(p, q, v) = - / u ■ mrdr = - ^ (piPj + ViV^G{q h qj) . (33) 

J i,j=l 

Substitution of the solution ansatz (|5Uj) for the momentum and its corre- 
sponding velocity ([3*2^1 into the radial equation © gives the system, 



E 



P. 



S(r - qi) . S'(r- qi) 



i.j 



5{r - gj) 5'{r - qj) 



5(r-qj) dG 
+ 2piPj -fr(r,qj) 



S(t — Qi^j Sir — Qi) dG 
+ (ViV 3 - ViVj) ^^G(r, qj) + v iVj ^^^r( r ' 9j) } = • 



Multiplying this system by the smooth test function rip(r) and integrating 
with respect to r yields dynamical equations for Pi and In particular, the 
ip(qi) terms yield 



Pi = -J2 ( Pi ^ff +ViVj { 



(dG{q u q j ) G(q u qj) 



dqi 



(34) 
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and, after integrating by parts, the ip'(qi) terms yield 



(35) 



By equation (j22J) we see that q^t) = f • u(gj,t), so the radius of the i— th 
cylinder moves with the radial velocity of the flow. 

This procedure is repeated for the component of the EP equation, by 
substituting the solution ansatz ()30I32|) into equation Q, to find the system 

/. 5(r- qt) . 5'{r - ft) 
> J Vi Viqt 



i.j 



8{r-qi) 8'{r-qi) 



5{r- qi ) fdG n G(r,qj) 
v iPj ( ^(r, Qj ) + 2^— 



0. 



Upon multiplying this system by rip{r) and integrating with respect to r, 
the term proportional to ^'(qi) again recovers exactly the g^— equation (|33j). 
The term proportional to ^(qi) gives 



(36) 



after using the q^— equation (j33|) in the last step. This integrates to 

Viqi = Mi = const, (37) 

where Mj are N integration constants. From the Hamiltonian viewpoint, this 
was expected: The angular momentum Mj is conserved for each individual 
circular peakon, because each circle may be rotated independently without 
changing the energy. 

Equations (J33j) and (f33j) for pi and qi may now be recognized as Hamilton's 
canonical equations with Hamiltonian, 

H(p, q, M) = ~ /u • mrdr = \ ^ (p iPj + MiMl^Q^ Qj ) . (38) 
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This is the same Hamiltonian as obtained from substituting the momentum 
map (|2"7|) into the kinetic energy KE in equation (jHSJ)- Hence, we may 
recover the reduced equations (JMj) for p i: f|H3|) for q i and (|37jl for M t from 
the Hamiltonian (J38j) and the canonical equations, 



dH . dH ■ dH n 

This result proves the following: 

Proposition 4.1 TTie momentum map \2T\) is a valid ansatz for rotating 
peakon solutions of the Euler-Poincare equation JTJj for geodesic motion. 



4.2 Solution properties 

The remaining canonical equation for the i— th Lagrangian angular frequency 
is, 

f)M M 1 

i< = m = ^m G(g " g ' ) = 7, !><')«<*.*)■ («> 

Thus, as expected, the ignorable canonical angle variables = with 
i = 1, 2, . . . , N, decouple from the other Hamiltonian equations. In addition, 
we see that 

q4i{t) = • u(q h t) , (41) 

so the angular velocity of the z— th cylinder also matches the angular velocity 
of the flow. Therefore, we have shown: 

Proposition 4.2 The canonical Hamiltonian parameters in the momentum 
map and solution ansatz \2T\ j provide a Lagrangian description in cylindri- 
cal symmetry of the flow governed by the Eulerian EP equation for geodesic 
motion (QJ). 



Angular momentum, fluid circulation and collapse to the center. 

Finally, the fluid circulation of the i— th concentric circle which is traveling 
with velocity u, may be computed from equations ()29|) and (J30|) (with a slight 
abuse of notation) as, 

/* f 27rA^ 
® m • dx = (p rm^dcf) = 2iiVi = -. (42) 

Jci(u) Jci(u) Qi 
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We see that the "angular velocity" Vi = Mi/qi is the fluid circulation of the 
i— th concentric circle. Since the angular momentum Mj of the i— th circle is 
conserved, its circulation Vi(t) varies inversely with its radius. Consequently, 
this circulation would diverge, if the i— th circle were to collapse to the center 
with nonzero angular momentum. 



5 Numerical Results for Radial Peakons 
5.1 Radial peakon collisions 

We consider purely radial solutions of equation (|SJ), with = 0, which 
satisfy 

dm 

- r = d r (rm r u r ) — m r d r u r . (43) 

at r \ / 

Such radial solutions have no azimuthal velocity. Without azimuthal velocity, 
the vector peakon solution ansatz (JSUj) for momentum reduces to the scalar 
relation, 

mr( r,t) = ± P ,(t) S ±^l. (44) 

i=l 

The corresponding radial velocity is 

N 

Ur(r,t) = ^2 Pi (t)G(r,qi(t)), (45) 

i=l 

where the Green's function G(r,qi(t)) for the radial Helmholtz operator is 
given by formula (|13|). Radial peakons of this form turn out to be the building 
blocks for the solution of any radially symmetric initial value problem. We 
have found numerically that the initial value problem for equation (|44j) with 
any initially confined radial distribution of velocity quickly splits up into 
radial peakons. This behavior is illustrated in Fig- 12 The initial distribution 
of velocity splits almost immediately into a train of radial peakons arranged 
by height, or, equivalently, speed. 

The head-on "peakon- antipeakon" collisions are of special interest. In 
the case of equal strength radial peakon-antipeakon collisions, the solution 
appears to develop infinite slope in finite time, see Fig. El This behavior is 
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Peakon iorma:ion ; rom inrial no n- rowing Gaussian 




Figure 2: The initial value problem: a Gaussian profile splits into peakons. 




Figure 3: Peakon-antipeakon collision of equal strength. 
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5 10 15 20 



Figure 4: Peakon-antipeakon collision of unequal strength. The smaller 
peakon's trajectory undergoes a large phase shift. 

also known to occur for peakon-antipeakon collisions on the real line. If the 
strengths of the peakon and antipeakon are not equal, then the larger one of 
them seems to 'plow' right through the smaller one. This is shown in Fig. 0] 

The figures shown were produced from numerical simulations of the Eule- 
rian PDE (J43|) . The momentum m r was advanced in time using a fourth-order 
Runge-Kutta method. The time step was chosen to ensure the Hamiltonian 
| J" m ■ urdr was conserved to within 0.1% of its initial value. The spatial 
discretizations ranged from dr = 10 -4 to dr = 0.02 depending on the desired 
resolution and the length of the spatial domain, and the spatial derivatives 
were calculated using finite differences. Fourth- and fifth-order centered dif- 
ferencing schemes were used for the first and second derivatives, respectively. 
The momentum m r was found from the velocity u T using the finite difference 
form of the radial Helmholtz operator, and the velocity u r was found from 
the momentum m r by inverting the radial Helmholtz matrix. For the peakon 
interaction simulations, the initial conditions were given by a sum of peakons 
of the form (|45jl for some chosen initial pi and qi. For a peakon collapsing to 
the center, which will be described next, the boundary condition at the ori- 
gin is important. If the PDE (j4Hj) were extended to r < 0, then the velocity 
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would be an odd function about the origin. In addition, when the peakon 
was sufficiently close to the origin (qi < 0.1), the sign of its momentum m r 
was reversed to begin its expansion away from the origin. 

For comparison with the simulations of the Eulerian PDE simu- 
lations of the Lagrangian ODEs were also performed. A fourth-order 
Runge-Kutta method was used to advance the system in time, and a time 
step was chosen to ensure the Hamiltonian (J38|) was conserved to within 
0.1%. The results of these simulations agreed with those of the Eulerian 
PDE simulations to within 1%. 

5.2 Bouncing off the center 

Let us first consider the case when only one peakon collapses onto the center 
with the angular momentum being zero. The Hamiltonian in this case is, 

H=^p 2 I 1 (q)K 1 (q), 
which can be approximated when q — > as 

H= {^ + o{q)^p 2 . 

Thus, the momentum p(t) is nearly constant just before the collapse time, 
t*, and is approximately equal to —2\frl, more precisely, 

p= -2VH + o(q). 

The equation of motion for q(t) yields, q = pIi{q)K\{q) = —\f~H + o(q). If 
q — > at t — > £*, then we necessarily have 

q{t) =VH{U-t) + +o((U - t) 2 ) , 

near the time of collapse t t*. 

The case of iV radial peakons can be considered similarly. If only one 
peakon (let us say, number a) collapses into the center at time t*, so q a (t) — > 
as t — > t*, and the motion of the peakons away from the center is regular in 
some interval (t* — 5, t* + 5) (as will be the case unless a peakon-antipeakon 
collision occurs during this interval), then conservation of the Hamiltonian 
implies 

P 2 a G(q a , q a ) + 2p a A + B = 2H, (46) 
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where 



.4 



B 



PiPj G (9i,Qj) 

(«J)A 



Since q a = min(qi, . . . , q^), the Helmholtz Green's function in expression 
(fT3~|) implies that the quantities A and B are bounded at times close to t*, 
and that G(q a , q a ) is bounded as well. Consequently, equation (jlUj) implies 
that p a is also bounded at times close to t — t*. 

Numerical simulations confirm our predictions: At the moment of the 
impact at the center, the amplitude of the peakon remains bounded and 
approaches the value of —\rH « —2.23, as illustrated on Fig. El 

The slope of the solution at the origin has to diverge. This can be seen 
on the example of a single peakon as follows. Since 



Thus, if q(t) — > as t — > t*, the slope d r u{r = 0,t) must diverge. Therefore, 
the following Proposition is true: 

Proposition 5.1 A radially symmetric peakon with no angular momentum, 
collapsing to the center, has bounded momentum and unbounded slope at the 
origin close to the moment of collapse. 

6 Numerical Results for Rotating Peakon Cir- 
cles 

Simulations of rotating peakons were performed for the Eulerian system of 
PDEs (JHEj) using the same numerical methods as those used for non-rotating 
peakons. Fig. El shows the results of an initial value problem simulation 
when u T is initially and is initially a Gaussian function. Radial velocity 
in both directions is almost immediately generated, and rotating peakons 
soon emerge moving both inward and outward but all rotating in the same 
direction. A rotating peakon approaches the center but turns around before 
reaching the origin. 



Or 



du 



p^I'MKM 



l + o(q)p(t) 
2 q(t) 
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Figure 6: An initial angular velocity distribution (with zero initial radial ve- 
locity) breaks up into rotating peakons which move both inward and outward. 
The radial velocity u r , the angular velocity u^, and the velocity magnitude 
lul are all shown. 
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Fig. shows a rotating peakon as it approaches the origin. A sort of 
angular momentum barrier is reached, and the peakon turns around and 
moves away from the origin. Thus, a peakon's behavior as it approaches the 
origin is reminiscent of Sundman's Theorem; if a peakon has nonzero angular 
momentum, then a full collapse to the origin will not occur. This result can 
be understood as follows. For a single rotating peakon, the Hamiltonian (J38|) 
becomes 

H= (V + ^J G(q,q) 

From the theory of Bessel functions we know that G(q, q) — * 1/2 when q — >• 0, 
and G(q, q) > for all q. Moreover, it can be shown that G(q, q) is strictly 
decreasing with increasing q > 0, so we conclude that G(q, q) > G(qo, go) if 
< q < q . Thus, when < q < q , 

M 2 

H>—G(q ,q ) 



q 2 



so 



g>max(g ,M^%Mj (47) 

The estimate (|47jl provides the lower bound q(t) can reach in the process of 
evolution for each value of the parameter go- This estimate can be further 
optimized as follows. Since G(q , g ) is strictly decreasing starting with the 
value of G(0, 0) = 1/2, there is q*(M, H) such that 



q , = M^/fe) (48) 

for each value of M, H > 0. Then, q > q*(M,H) is the desired optimal 
estimate. We can summarize this result in 

Proposition 6.1 Consider a rotating radial peakon with given values M, H > 
0. If we define q*(M, H) by fifify , then in the process of evolution this peakon 
cannot be closer than q*(M,H) to the origin. 

Note One can show that g*(M, H) ~ M/V2H + O (M 3 ) for M -»■ 0. 
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Figure 7: A rotating peakon turns around near the origin. The radial velocity 
changes sign as the angular velocity reaches a maximum. 
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7 Conclusions 

The momentum map ()16|) for the action of diffeomorphisms on (closed) curves 
in the plane was used to generate invariant manifolds of measure- valued solu- 
tions of the Euler-Poincare equation (0) for geodesic motion on the diffeomor- 
phism group. The H 1 norm of the fluid velocity was chosen for the kinetic 
energy of a class of solutions whose momentum support set was a family of 
N closed curves arranged as concentric circles on the plane. These finite- 
dimensional invariant manifolds of cylindrically symmetric solutions general- 
ized the N— peakon (soliton) solutions of the Camassa-Holm equation from 
motion of points along the real line, to motion with cylindrical symmetry. 
This cylindrical motion included rotation, or, equivalently, circulation, which 
drives the radial motion. The momentum map with non-zero circulation for 
these concentric circles yielded a generalization of the circular CH peakons 
that included their rotational degrees of freedom. The canonical Hamiltonian 
parameters in the momentum map and solution ansatz (|27|) for the concentric 
rotating circular peakons provided a finite-dimensional Lagrangian descrip- 
tion in cylindrical symmetry of the flow governed by the Eulerian EP partial 
differential equation for geodesic motion (0). Numerically, we studied the 
basic interactions of these circular peakons amongst themselves, by collisions 
and by collapse to the center, with and without rotation. 
The main conclusions from our numerical study were: 

• Collapse to the center without rotation occurs with bounded canonical 
momentum and with vertical radial slope in velocity at r = 0, at the 
instant of collapse. 

• For nonzero rotation, collapse to the center cannot occur and the slope 
at r = never becomes infinite. 

The main questions that remain are: 

• Numerical simulations show that near vertical or vertical slope occurs 
at head-on collision between two peakons of nearly equal height. A 
rigorous proof of this fact is still missing. 

• Is the motion integrable on our 2N dimensional Hamiltonian manifold 
of concentric rotating circular peakons for any N > 1 and any choice 
of Green's function? 
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• How does one determine the number and speeds of the rotating circular 
peakons that emerge from a given initial condition? 

• How does the momentum map with internal degrees of freedom gener- 
alize to n dimensions? 

All of these challenging problems are beyond the scope of the present paper 
and will be the subjects of future work. 
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8 Appendix: Extension of one-dimensional 
peakons 

In this appendix, we show how to obtain the momentum line peakons, which 
are generalizaion of the line peakons. The standard ansatz for the regular 
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:,t) = ^p i (t)6{x-q i (t)), (49) 



one-dimensional peakon is 

N 

m{x, '] 

8=1 

where m satisfies the one-dimensional version of (JTJ). We propose the follow- 
ing extension of these solutions: 

N 
i=l 

where x, y are unit vectors in x, y directions, respectively. The solution lives 
on line filaments, which are parallel to the y axis and propagate by translation 
along the x axis. However, the y component of momentum now has a non- 
trivial value. Such solutions represent momentum lines which propagate 
perpendicular to the shock's front and "slide" parallel to the front, moving 
surrounding 'fluid' with it. Upon substituting (|50|) into equations of motion 
(JTJ), we see that the x and the y components of (HJ) both give the same 
equation of motion for qi(t): 

3 

This compatibilty is what makes the factorized solution (150(1 possible. The 
equation of motion for pi is 

PM = ^2(PiPj + ViVj)G'(qi, qj) , 
j 

and for v i: 

i>i = 0. 

Thus, Vi can be considered as a set of parameters. The (p,j, qi) still satisfy 
Hamilton's canonical equations, with Hamiltonian now given by, 

1 N 

H = - (PiPl + V i V i) 1j) ■ ( 51 ) 

Finally, the "angle" variables yi(t) conjugate to Vi{t) with canonical Poisson 
bracket {yi , Vj} = 5ij satisfy: 



